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Abstract 

Wc prove that, in a generic single- field model, the consistency relation for the 3-point function in the squeezed 



limit receives corrections that vanish quadratically in the ratio of the momenta, i.e. as (fci/feg)^. This implies 

Q I that a detection of a bispectrum signal going as l/k\ in the squeezed limit, that is suppressed only by one 
5— ( ' 

-4-^ ' power of k^ compared with the local shape, would rule out all single-field models. The absence of this kind of 
C/3 ■ , 

' terms in the bispectrum holds also for multifield models, but only if all the fields have a mass much smaller 
than H . The detection of any scale dependence of the bias, for scales much larger than the size of the haloes, 

' would disprove all single-field models. We comment on the regime of squeezing that can be probed by realistic 

CSl ' surveys. 



o 



1 Introduction and main results 



The understanding of a physical system is often enlightened by taking parametric limits. In the case of the 
primordial 3-point function of scalar perturbations, it is useful to consider the so-called squeezed limit, when 
one of the momenta in Fourier space becomes much smaller than the other two. From the theoretical point 
of view, it is well known that in this limit a consistency relation between the bispectrum and the power 
spectrum can be written, with the only assumption that perturbations are generated by a single dynamical 
degree of freedom [1, 2, 3]. More recently the squeezed limit has become of experimental interest, since one of 
the most promising signature of primordial non Gaussianity in the large scale structure, the scale-dependent 
bias [4, 5, 6], is only sensitive to the squeezed limit of the 3-point function. Motivated by this experimental 
interest, in this paper we study the corrections to the consistency relation as we slightly depart from the 
squeezed limit. 

The consistency relation implies, in practice, the absence of terms diverging in the squeezed limit like 
fc^'^, where k^ is the small momentum that we send to zero in the squeezed limit. Single- field models on the 
market typically give a bispectrum suppressed in the squeezed limit by two powers of fc^, i.e. diverging only 
like l/kh- A natural question is therefore whether this is a general result, or also an intermediate behaviour 
like \/k\ is possible in single-field models. This question is particularly relevant for measurements of the 
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scale-dependent bias; while a behaviour l/fc^ is impossible to observe, being completely degenerate with the 
standard bias, a l/fc£ behaviour would leave interesting and potentially distinguishable signatures. 

In this paper we show that it is impossible, in any single-field model, to have a l/fc|_ divergence of the 
bispectrum in the squeezed limit. More precisely we will modestly sharpen up the consistency relation, 
rewriting it as 



the short modes ^2 and fc3, there are no linear corrections to the consistency relation. We will show this in 
Section 2, reviewing first of all the general derivation of the consistency relation, that directly leads to Eq. (1) 
and then, in the remaining of the Section, checking that all corrections are at least quadratic in the long 
wavelength momentum fci . 

Our results arc valid for any single-field model, with the only requirement that the unperturbed history is 
an attractor. No slow-roll approximation is required and indeed in Section 3 we verify explicitly the absence of 
linear corrections to the consistency relation in models with features or modulations in the inflaton potential. 
It holds also in models where scale invariance of the spectrum is not associated to a quasi de-Sitter phase, 
like models with a strongly varying speed of sound [7, 8] or for the so-called adiabatic ekpyrosis [9, 10], as we 
discuss in Section 3. For simplicity, in the following, wc will always call the scalar field infiaton, although our 
conclusions arc not restricted to an inflationary scenario. 

It is rather intuitive why the consistency relation only receives quadratic corrections. In a single-field 
model a long wavelength mode docs not change the local physics ~ it is adiabatic - and in particular cannot 
change the short-scale 2-point function. This is true as long as we can neglect the curvature induced by this 
long mode, i.e. up to order fc£. Of course what we mean by 'cannot change the short scale 2-point function' 
depends on the gauge we choose (see for example [11]): indeed the 3-point function of C is not vanishing in the 
squeezed limit, but it is proportional to the tilt of the spectrum, since the background mode induces a (locally 
unobservable) rescaling of the coordinates. In this paper we show the linear corrections are completely absent 
if the consistency relation is written in the form (1). 

Multi-field models do not satisfy any consistency relation. Indeed it is typical of multi-field models to 
have a local divergence of the 3-point fimction in the squeezed limit, i.e. l/k\, with a generic, potentially large 
parameter /^l^' replacing the small tilt in eq. (1). However, in Section 4, we argue that also in multifield 
models it is quite difficult to have linear corrections to the local form, i.e. going as l/fc|. In particular we 
will show that the absence of this kind of term is general if all the fields participating in the generation of 
primordial perturbations are nearly massless, i.e. m <^ H. No general conclusion can be drawn for models with 
fields with masses of order Hubble during inflation, as for example in the so-called quasi-single field models 
[12]. As the consistency relation applies for k^/ks — > 0, in Section 5 we quantify the regime of squeezing that 
can be realistically probed by measurements of the scale dependent bias. Conclusions arc drawn in Section 6. 

2 Single-field consistency relation with subleading corrections 

In this Section we review the standard arguments leading to the consistency relation, discussing the necessary 
assumptions and approximations, and show how these can be partially relaxed. In the standard derivation 



(C(fci)C(fc2)C(fc3)) ^ -(27r)3<5(fci +k2 + k3)P{ki)Piks) 



dln(fc|P(fcs)) 
din ks 




ki < ks , (1) 



where ks = {k2 — fc3)/2. Once the bispectrum is written in terms of fcg, i.e. symmetrically with respect to 
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[1, 2, 3], the long-wavelength field is assumed to be a spatially homogeneous background, and the resulting 
statement is therefore valid strictly when fci — >■ and k2 — > —k^. We will show that one can include mild 
spatial variations of the background, without substantially modifying the result. The main result of this 
Section, and of the paper, will be to prove that the bispectrum written in terms of ks = (^2 — k3)/2 does not 
receive linear corrections in ki/ks when moving away from the exact squeezed limit. 



2.1 Review of the standard leading order computation 

We review here the proof of the leading order consistency relation, following [3]. In the squeezed limit, the 
mode C(fci) leaves the Hubble radius much earlier than the other two, and it acts as a constant background 
for the other two modes. The bispectrum can be computed in two steps: first one computes the two-point 
function (C2C3)f(a;) in the presence of a background C,{x), and then averages this two-point function over the 
realizations of the background: 



(a(fci)c(fc2)c(fc3)) ^ (C(fci)(c(fc2)c(fc3)). 



(2) 



We therefore need the function (C(^2)C(^3))f(a;)j which can be conveniently computed in coordinate space 
before taking a Fourier transform. When the background is assumed to be homogeneous, C{x) = C it can 
be reabsorbed simply by rescaling i' = e^a;*. After such rescaling C, no longer appears in the action, so that 
the 2-point function in the new coordinates is the same as in the absence of C,. In the limit of ^ constant, 
this reads {C,{x2)C,{xz)) q = (C(x2)C(2^3))- When C, is slowly-varying, it can be conveniently evaluated in the 
middle point = {x2 + a;3)/2 to get x^ ~ X2 — x^ ~ X2 C(^+)(2?3 ^ ^2)- Therefore one can write the 2-point 
function at linear order in Q as 



(C(x2)C(^3)>lc» - ^(^3 - .T2) +C(^+)[(^3 - X2) ■ V^(|^3 - X2\ 



(3) 



where £,{\xz — a?2|) = {2tt)^^ j 6?kse'^''^'^^'-^~^'^^ P{ks). It is important to notice that corrections to this ex- 
pression do not contain terms proportional to V^, since there is no vector to contract with: everything is 
symmetric under the exchange X2 <^ 2^3, so we cannot define a preferred direction. It is easy to verify that 
corrections of order d^C, will correct the final result only at 0{k^/kg). The 3-point correlation function will 
then read 



(C(2?l)C(^2)C(^3)) (C(.fi)C(f+))[(x3 - X2) ■ VC(|f3 - £2\)] 

d^kL f d^ks 



(27r)3 J (27r) 



P{kL)P{ks) 



dks 



(4) 



where X- = .T3 — X2. After integrating by parts, inserting 1 = J d^ki S{ki + k^) and using the relation 
{d/dks) ■ [ksP{ks)] - P{ks)d\n[klP{ks)]/ dinks one gets 



d^fci d^fcL d^fcg 
(2^0^ 



i27rfd{k,+kL)P{ki)P{ks) 



dlnfc|P(fcs) 



d In ks 



(5) 



In this expression, one can set /c^ = A:2 + fc3 and ks = {k2—k^)/2 and get — ifc/^-.T+-|-iA:5-x_ 
Changing variables in the integration and Fourier transforming, one finally obtains 

din k%P{ks) 



(a(fci)C(fc2)C(fc3)> - -(27r)-^<5(fci 



■ k3)Pik^)Piks)- 



din ks 



-ik2-X2-ik3-X3. 



(6) 



which is the leading order result of (2). 



3 



Notice that writing ks in terms of ^2 and it is straightforward to get P{ks) = ^[P{k2) + Piks) + 
0{kL/ks)'^], where the hnear corrections cancel out. Therefore, since the logarithmic derivative above is just 
the definition of the tilt of the power spectrum, the consistency relation can also be written as 

(C(fci)C(fe)C(fc3)> = -i^TrfSih + fc2 + h) {ns ~ 1) P(fci) ^(^^^+^(^^) (7) 

up to corrections that are quadratic in k^/ks- In this way it is explicit that the behaviour in the squeezed 
limit is the same at the one of a local shape up to quadratic corrections. 



2.2 A non-standard proof with standard methods 



The fact that the effect of a background mode ( amounts to a rescaling of the coordinates can also be seen at 
the level of the action. In the exact squeezed limit, the long-wavelength field is constant in time and space: 
the only relevant cubic terms are thus the ones with (at least) one leg without derivatives, which will be 
evaluated on the long-wavelength mode. If we focus on these terms in Maldacena's cubic action for ( - see 
eq. (3.7) of [1] ~ after some integrations by parts one ends up with 



'S'2 + S3 



Ml, 



A^xea^ (1 + 3C)C'-(1 + 0(9,0 Va' 



(8) 



where C, is the leg that will be evaluated on the long mode and e = —H/H^. In this expression we have also 
written the quadratic action for C, as it makes evident how the consistency relation works: Q can be got rid of 
through a spatial coordinates rescaling — )- a;'(l -I- and we are left with the unperturbed quadratic action, 
so that the 2-point function calculation is insensitive to the background wave in the new coordinates. Exactly 
the same will happen for any single field model: the same coordinate change eliminates the cubic interactions 
in the squeezed limit. 

Although the rescaling argument is very intuitive, for completeness we give an alternative proof using the 
standard computation of the 3-point function from the cubic action. For simplicity in the remainder of the 
Section we focus on the simplest slow-roll inflationary scenario. For a different derivation of a very similar 
result see [13, 14]. Using the in-in formalism, the tree level bispectrum reads 



(C1C2C3) = I 



C1C2C3, / d^xea^C (3C'-(V0Va') 



(9) 



where the coordinate dependence of the fields Q = C{t, Xi) and C = 1 2^) is understood. Under the assump- 
tions above, we can retain only the contractions of Ci (the long-wavelength mode) with C, and neglect those 
with C and V(^. Moreover, in the limit in which we are working, (CiC) is just a real constant. We therefore 
have 



(C1C2C3) / d3.T(CiO( 



C2C3 



At' e a' 



(3C^ ~ (VOV«') 



(10) 



We now notice that if e were not there, the integrand would be exactly the time derivative of the Hamil- 
tonian of a free scalar field in de Sitter. We try therefore to express the integrand in terms of the quadratic 
Hamiltonian and Lagrangian densities '^2 = a'^e (^'^ + (VC)^/a^^ and £2 = a'^e (^t^ — (VC)^/a^^ , in the form 
\f{t)T-L2] + g{t)C2- The reason is that, once this is done, the integral of the first term simply gives f{t)H2. 
Since H2 generates the time evolution of the free fields, its commutator with C2C3 is proportional to (C2C3) 
and therefore vanishes outside the Hubble radius. We will therefore only have to compute the integral of the 
second term, 6t'g{t')C2- 
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The integrand can be expressed in terms of 'H2 , ^2 and their derivatives as 



ah 



3C^-(VC)V«' 



■H2 + 2C2 



3C - (VC)Va^ 



1 

H 



•H2 + -C2 
e 



Writing '^2 /H ~ {H2 / H) — e'H2 and combining with the two previous expressions one gets 



ah 



3c' - (ycf/a 



2 /„2 



1 



H 



He 



2e £2 



(11) 



(12) 



1 - € 

The same process can be repeated to bring 1/(1 — e) inside the total time derivative (as done with 1/ H), and 
iterated up to arbitrary order, so that in general one has 

1 



a\ 



3C^-(vc)7«' 



2e + 0(e") C 



(13) 



Xl-e)H 

As explained above, the commutator with C2C3 of the integral of the first term vanishes outside the Hubble 
radius, and we just have to integrate the second (which also vanishes in the scale invariant limit). In order 
to do this, it is convenient to express £2 as £2 = (a'^eCC) ~ ^ ■ (c'eC^Oj using the equation of motion for the 
free field. The second term can be integrated by parts to give V(C(a;i)C(x)), which can be neglected under 
the above assumptions. All we have to compute is thus 

At' (-^+2e + Oie')] (ahCC) =(^+2e + 0{e^)) ahCC- f di'a^e f + 26 + 0{e^)] CC , (W) 



He 



He 



He 



where the second integral is already of order O(e^) and can be neglected. Recalling that the conjugate 
momentum of C is tt = 2a^eC, the three-point function thus reads 



(C1C2C3 



2e + 0{e^)^ j d3x(CiC(a'))([C2C3,^(^)C(a')]), 



2 \He 

where the fields are now all evaluated at the same time. Using the canonical commutation relations [Ci, C(-'^^) 
and [C,iT'K{x)\ = iS{xi — x) one finally has 

(C1C2) + (C1C3) 



(15) 



(C1C2C3 



He 



+ 2e + 0{e') 



-(C2C3 



(16) 



The proof is then completed by showing that e/{He) = 2(e — S), with S = (p/{H(j)), so that e/{He) + 2e = 
4e- 2(5 = -{ns - 1). 



2.3 Including gradients: proof of the absence of linear corrections 

We want to prove that, in any model, there are no linear corrections to the consistency relation, eq. (1). 
The derivations of the consistency relation above consider only interactions in the cubic action for C that 
do not have derivatives on the background C and, moreover, take ^ as time-independent. Given that in 
this approximation we get eq. (1), to conclude we need to show that all terms that we neglected give only 
corrections of order ©(/ci/fcs)^ or higher. In doing so we have also to consider what comes from the solution 
of the constraints which, introducing non-local terms, could invalidate our arguments. 
To summarize, we need to prove the following statements: 

1. Corrections will arise when we take into account the time-dependence of the background mode, which 
above has been taken as constant. These subleading contributions are 0{k^) as we will see in Section 
2.4. 
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2. The cubic action will contain terms with a single spatial gradient acting on the background mode, di(^, 
which naively give 0{kL) corrections. However this behaviour cancels after symmetrization: the vector 
fci, from the derivative will be contracted in turn with k2 and fcs, which are equal and opposite at 
leading order in the squeezed limit. In other words, the short-mode 2-point function has no direction - 
as it is symmetric if we exchange the 2 points - so that the background mode gradient has nothing to 
contract with^. Thus corrections are 0{k\). 

3. Terms in the action with time derivatives on the long-wavelength mode, C, give 0{k^) corrections. This 
is a consequence of point 1, whenever ( tends to a constant. 

4. Some cubic interactions originate from solving SN and Ni through the (linear) constraint equations. 
In Section 2.5 we prove that 6N = 0{k\), and that Ni is 0{kL) but, following point 2 above, it gives 
0{k\) corrections after symmetrizing over the short modes. 



2.4 The decay of the decaying mode 

At lowest order in derivatives, the most generic single-field quadratic action for C is 

(17) 



where we have allowed an arbitrarily time-dependent speed of sound Cg. After defining a modified conformal 
time dy = {cs/a)dt, this action reads [8] 



M^iJ dyd^xq^ 



(18) 



where primes stand for derivatives with respect to y and = a^ejcs- The equation of motion from this action 
is ^ 

C" + 2^C'- v^C = o. (19) 

In the long wavelength limit, i.e. up to terms quadratic in the gradients, this equation admits two solutions, 
Cfe ~ constant and C,k ^ J dy{l/q^). Assuming that the constant solution is dominant - which is equivalent to 
requiring that the background solution is an attractor - it is easy to realize that (approximate) scale invariance 
implies (approximately) q oal/y [8, 15]^. This means that for the second solution ocy^: the decaying mode 
dies off as y'^ in all models which give rise to an approximately scale-invariant spectrum. This behaviour 
may be changed locally if we allow for a feature in the power spectrum, or being only valid on average when 
a periodic modulation of the power spectrum is considered. In the squeezed limit of the bispectrum, the 
long-wavelength mode freezes {ky ~ 1) much before the others and the decaying mode dies off as y^ after 
that. Therefore, when it is the turn of the short modes to freeze its amplitude is suppressed by (k^/ks)^, 
much smaller than the corrections we are interested in. In particular, as we can neglect the decaying mode, 
^ ~ 0{k\) as it comes from the gradient corrections to the dominant mode. Therefore also C, terms acting on 
the long- wavelength mode cannot give any correction linear in kj^/kg. 



^This argument does not work for the 3-point function of tensor modes, as in this case the polarization tensors 
introduce additional structure. Indeed one can check that there are linear, 0{kL), corrections to the consistency 
relation for tensor modes discussed in [1]. 

^We will come back to discuss various ways to achieve this in Section 3.2. 
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As a possible caveat to the logic above, one could imagine that the decaying mode, which decays after 
freezing, is boosted by some effect inside the horizon which makes the mode depart from the Bunch-Davics 
vacuum. This happens for example in the presence of sharp features of the potential. However it is easy to 
realize that this cannot change our conclusion, at least for sufficiently squeezed triangles. As we squeeze the 
triangle the out-of-the-horizon suppression (/c^ /ks)^ becomes smaller and smaller. For the decaying mode to 
still give some effect, and in particular to give linear corrections to the consistency relation, we would need 
to boost without bound its amplitude before freezing. This is impossible, as at a certain point the linear 
perturbation theory will break down. 

The quadratic action (17) is only valid for models without higher derivatives. However, there are inter- 
esting models of inflation - ghost inflation and its generalizations [16, 17] - for which the spatial kinetic term 
has four derivatives: (dfC,)'^. Following the logic above, with an appropriate redefinition of the time variable, 
one can cast the quadratic action in the form 



(20) 



The y dependence in front of the spatial kinetic term is inspired by the case of Ghost Inflation: in a de 
Sitter background, the additional derivatives of the spatial kinetic term give a factor 1/a^ — t'^ compared to 
the time kinetic term. Written in this way, it is easy to realize that scale invariance (and the requirement 
that ( goes to a constant out of the Hubble radius) is achieved with z oc 1/y. Indeed in this case one has 
a "fictitious" de Sitter dilation symmetry (x, y) — >■ X{x,y), which implies the scale invariance of the 2-point 
function. Given that the discussion above does not depend on the spatial kinetic term and we have the same 
time dependence in front of C'^ , the conclusion about the behaviour of the decaying mode outside the horizon 
remains unaltered. 



2.5 Constraints on the constraints 



Let us look at the solution of the constraint equations in general models. To compute the bispectrum, we just 
need the solution for the lapse SN and the shift iV' at linear order [1]. 
Consider the general action for inflation 

Mi 



5, 



EH 



^Pl 



Vh 



Sn 



(21) 



where Eij = \{hij — ViNj — VjiV^) and Sm is the inflaton action. 

We want to understand the general structure of the linearized constraint equations: 



PI 



i?(3) _ N-'^{EiEf - E^ 



_5S__^/h 

0^^^VhA4,V,[N~\Ei-SfE) 



SN 



(22) 
(23) 



For single-field models, the most general expression for the inflationary Lagrangian is given by the effective 
theory of inflation [18]. In the comoving (or unitary) gauge, the second order Lagrangian reads: 



M^.Hg"'' - MUZH^ + H) + iM4(500 + if 



00 



l)SE - ^APSE^ - ^AIHE]5EI 



(24) 
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where the mass scale M and the dimensionless coefficients Ci are in general time-dependent and the dots stand 
for higher derivative terms we will discuss later. 

In this gauge, where the scalar field is unperturbed and hij ~ a^{i) exp (2C) 5ij, we have at linear order 

i?(3) = , (25) 

5E] = C-^;. - [d'Nj + djN')/2 , (26) 

6E = 3C- diN' . (27) 

The constraints therefore read 

25N{iMl^H'^ + Ml^H - 2M^) + ^^N'' {2MI^H + cjA/^) = -2MI^—^ + ( {QM^^H + SciA/^) , (28) 

a,<5iV(2i7 A/|i + ci Af 3) + d^djN^ (c2 + ca)^/^ = a,C(2M|i + SczAf ^ + cg A/^) , (29) 

where we have used that the shift must be a gradient of a scalar function, = diip, as there are no vector 
modes. Eq.s (28) and (29) are inhomogeneous linear equations for 6N and diN^ sourced by terms proportional 
to C and As we discussed in the previous Section, in the long wavelength limit ( ~ 0{kj^). The equations 
thus imply that both 5N and diN"^ are 0{kj^) in the long wavelength limit. This means that neither will give 
corrections 0{k]^) to the squeezed limit of the 3-point function: 6N can be neglected, while for TV* - which is 
0{ki,) - one can apply the same argument we used for diC, above: linear corrections cancel once symmetrized 
over the short modes. 

What about higher derivative terms left out from the action (24)? These can be built adding suitably 
contracted spatial derivatives V,; or derivatives orthogonal to the time surface g'^^V^ (in this case the upper 
index zero can remain unpaired.) For example let us look at the operators 

y (V,r)' + y (5°^V,r)' + d.g^^V.g^'SE , (30) 

which will generically appear when considering inflaton Lagrangians with two (or more) derivatives acting on 
the inflaton £((/>, V^V^^), as for example in [19, 20]. The first operator just contributes to the linearized 
Hamiltonian constraint with a term —AdiV^SN, which is suppressed by fc£ with respect to the others. In 
general additional spatial derivatives just contribute to corrections suppressed by additional powers of k\ to 
the solution of the constraint equations and thus do not change our conclusions. The second operator is more 
interesting because it contributes to the Hamiltonian constraint with a term —4(^2 + 3Hd2)SN — Ad25N . The 
Hamiltonian constraint becomes a dynamical equation for 5N . Of course, these higher derivative corrections 
must be treated perturbatively, as in any effective field theory, substituting into them the solutions of the 
constraints at lowest order. In this way, they cannot change the fact that in the long wavelength limit both 
5N and diN'^ are 0(k\). The same arguments applies to the third operator proportional to d-^ and in general 
to all higher derivative operators. 

3 Explicit check in various models 

It is straightforward to check cq. (1) in various models. We did it explicitly for the standard slow-roll result of 
Maldacena [1] and when one considers a small deviation from unity of the speed of sound [21]. It is well known 
that, neglecting slow-roll corrections, higher-derivative terms only give equilateral shapes, i.e. suppressed by 
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/c£ in the squeezed limit [22]. The absence of hnear terms should be valid also including slow-roll corrections, 
though the explicit check in the formulas of [23] is challenging. In the following two Sections we comment on 
two other classes of single- field models which satisfy our equation (1) in an interesting way. 



3.1 Modulations and features in the inflaton potential 

Departures from slow-roll represent, together with higher derivative terms, a possible source of large non- 
Gaussianity. These may appear as periodic modulations of the inflaton potential [24, 25] or as localized 
features [26, 27, 28]. 

Reference [25] gives an analytic expression for the bispectrum produced by models with a periodic mod- 
ulation of the inflaton potential 



(C(fcl)C(fc2)C(fc3)) (X 



l,2l,2yL2 



sm 



/2^ 



log kt/k^ 



/2i; 



log kt/h 



kt{kl 



kik2k3 



/2i; 



log kt/k^ 



(31) 



where kt = fci + ^2 + fca and fc» is a pivot scale. The result is obtained as an expansion in //-^/Se* ^ 1, i.e. in 
the limit of many oscillations per Hubble time and the dots stand for higher order terms in this expansion. 
Reference [25] shows that this expression satisfies the consistency relation at leading order. Notice that in the 
squeezed limit the second term in brackets dominates, though suppressed by f /\/2e^, as the first term diverges 
as 1/kf^. This behaviour of the first term, however, seems to blatantly violate our thesis. Terms subleading 
in //\/2e* luckily come to our rescue: they become relevant in the squeezed limit as the argument of the 
trigonometric functions contains \/2e^/ f. Indeed one can explicitely check that the second term cancels the 
l/fc£ divergence at zeroth order in //^/2e*, and the cancellation also works at first order in f/y/2e^, taking 
into account the third term in the expression. Though higher order terms have not explicitly been worked 
out, we expect this kind of cancellation to occur at any order. 

It is important to notice that eq. (1) applies only for "very" squeezed triangles 



kL 



/ 



/2i; 



< 1 



(32) 



and it is not enough to require fc^ <ti kg. This makes sense physically: the modes go through a resonance 
when they are a factor //■\/2e, shorter than the Hubble radius, i.e. when they have a frequency comparable 
to the one of the background. To deduce the consistency relation, we need the long-wavelength mode to act 
as a background for the others already at the time of the resonance. This is a stronger requirement than in 
slow-roll models, where it is enough that the long mode is frozen when the others cross the Hubble radius. 
When //v'2e* is very small, the consistency relation holds only for extremely squeezed triangles. In this case 
it is important to check how squeezed is squeezed for a particular observable, as the relevant range of scales 
may not be described by our eq. (1). We will address this point regarding the scale-dependent bias in Section 
5. 

Similar considerations apply to the case of localized features. In [28] an analytic expression for the 
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bispectrum is obtained using the generalized slow-roll approach which reads 



(C(A:i)C(fe)C(fc3)) « 



htO htO A^O 



with 



h = 
h = 



Ioih)kik2k3 + h{kt) kfkj + l2{kt)kt{kl + + k 



/ — G'B(logr)(A:tT)sin(A:tr) 
G'g{\ogT) cos{ktT) 
Gb (log r)(/ctT)— 

fc(T 



(33) 



f 


dr 


/ — oo 


T 


f 


dT 


J —oo 


T 


f 


dT 


/ — oo 


T 



(34) 



where Gs is a function of the background quantities and thus encodes the feature of the potential. Expanding 
kt around Iks inside the integrals, it is straightforward to check that terms going as l/fc|, cancel among the 
three contributions. As in the previous case, the expansion works only when k^ is sufficiently small that the 
mode is already frozen while the short ones are still in the vacuum and not yet perturbed by the features of 
the potential. 

Another situation in which the consistency relation only holds asymptotically, for very large squeezing, is 
with modified vacuum states; see the discussion in [29]. 



3.2 Non-standard single-field models with scale invariant spectrum 

In the presence of a time-dependent sound speed Cg [7, 8] or a rapidly varying equation of state [9, 10], one 
can obtain a scale-invariant spectrum even when the background evolution is very different from a quasi de 
Sitter expansion. 

In general, given an action of the form (18), the power spectrum will be scale invariant if and only if 
((' I q = To sec this, write down the equation of motion for the variable Vk = q^k 

d^vk + (^k^ - j^k ^ , (35) 

which will be the same as in dc Sitter when q" /q = with the only difference being that here the time 

variable is y. In general this can be realized if q cx l/y (case I) or q oc (case II). In case II the growing 
mode solution gives C (x vy~'^ which goes like y~^ for modes much larger than the sound horizon \ky\ <C 1: 
the background is thus not an attractor and it is quite contrived to get a viable scenario [15]. These models 
do not satisfy any consistency relation and we neglect them in the following. 

In case I, the linear equation of motion for C derived from the action (18) is exactly the same as in usual 
slow-roll inflation with the only difference being that the time variable is y, and therefore one expects C to 
have a similar behaviour: it is constant outside the sound horizon \ky\ <C 1 with corrections being at least 
quadratic in k^, and the decaying mode of this solution will be suppressed by y^, as explained in section 
2.4. One therefore expects all the arguments of previous sections to apply to these models, with one notable 
difference: scale invariance of the spectrum does not imply a scale invariant bispectrum. These models 
will exhibit in general a strongly scale dependent 3-point function, which implies that they become strongly 
coupled after a limited (~ 10) number of efolds [15]. Let us now focus on specific realizations of this case. 
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If one assumes that e is constant, an expanding background has q = —l/y, which for constant sound speed 
corresponds to de Sitter expansion air) = — l/r. However, a varying Cs allows for sizeable departures from de 
Sitter expansion [7, 8] . For instance, if e = const, a scale invariant spectrum is obtained for es = — 2e, where 
Cs = esHcg. Of course, the requirement of having inflation imposes e < 1. A stronger bound comes from 
the running of the gravitational wave spectrum, which does not depend on the speed of sound, and which 
translates into e < 0.3. Non Gaussianities can nevertheless be large and peaked in the equilateral limit since 
the speed of sound is small. Since the sound speed decreases with time, and will be smaller when smaller 
scales cross the horizon, the bispectrum will be larger at smaller scales. 

Similar considerations apply to the contracting "adiabatic ekpyrosis" scenario of refs. [9, 10]. In such a 
model ekpyrosis is driven by a single scalar field, and perturbations are generated during a phase in which 
the universe is slowly contracting, i.e. a{T) « constant, but the "slow-roll" parameter e varies rapidly with 
time, e ^ 1/t^. The quadratic action for ( is again given by an action of the form (18) with a speed of sound 
equal to unity and q (x 1/y. Due to the fact that e increases with time, non-Gaussianities in this model will 
be strongly scale-dependent and are expected to be larger at small scales. 

Even though the bispectrum is not scale- invariant, since we are dealing with single-field models the 
"squeezed limit" fci at fixed kg is expected to be of the form 



(C(fci)C(fc2)C(fc3)) = ~i27TfSiki +k2 + k:i)P{kL)P{ks) 



dln(fc|P(/cs)) 2, 



(36) 



d In ks 

Due to the absence of scale invariance in the bispectrum, we cannot say that the corrections to the squeezed 
limit will behave as 0{kj^/kg), but only that they will be proportional to at fixed ks- The scaling in ks 
will depend on the evolution of the speed of sound and "slow-roll" parameters with time, which fixes the scale 
dependence of non-Gaussianities. 

One can check explicitly that, in these models, the bispectrum does indeed take the form of equation (36). 
The case of varying speed of sound is discussed in [8] and the bispectrum is given in their eq. (7.38)^. Using 
those expressions, it can be easily verified that the corrections to the squeezed limit are quadratic in k^ for 
each operator of the action except for ('('^ and C(i9C)^- each one produces a bispectrum that has zeroth and 
linear order contributions in /cl, even for perfect scale invariance, but their contributions cancel. 

The non-Gaussianity in the adiabatic ekpyrotic model was computed in ref. [10]. In order to compute 
the bispectrum, the authors perform a perturbative expansion in a parameter 1/c <C 1, which must be small 
in order to guarantee a sufficiently wide spectrum of scale-invariant perturbations. One can check in their 
section 5.1 that the leading contribution to the bispectrum in the squeezed limit satisfies eq. (36). Subleading 
contributions, given in their section 8.1, are not of this form in the squeezed limit. However, since they are 
subdominant in 1/c one expects the full bispectrum at subleading order to receive contributions from terms 
ignored in the computation of the leading piece. The total bispectrum at each order in 1/c is expected to 
satisfy our results. 



4 Multi-field case 

How much of what we said above can be extended to multi-field models? Of course multi-field models do not 
satisfy a consistency relation for the 3-point function. Indeed, the SN formalism [30] relates, locally in space, 

''There is a typo in the published version of the paper of their expression A^Ji\ which has been corrected in the 
latest preprint version. 
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the fluctuations of tlie scalar perturbations in the spatially flat gauge to the final curvature perturbation 

C(f) = AiSifi'ix) + BijS^'{x)6ip^{x) + ... (37) 

This introduces some local non-Gaussianity for whose amplitude is no longer related to the tilt of the 
spectrum: the leading term in the squeezed limit is now a free parameter. However, one can still wonder 
whether it is possible to have linear corrections. In the following, we will show that linear corrections are 
absent only if one makes the additional assumption that all the fields have small mass (compared to H), and 
therefore a spectrum close to scale invariancc^. 

Before doing so, notice that, contrarily to the single-field case, here we have no control over the slow-roll 
corrections due to the deviation from de Sitter and the mass of the fields; we can only assume that they are 
negligible. When these corrections become large, i.e. in the presence of fields with mass of order H, we cannot 
say anything general. Indeed, in the presence of large deviations from scale-invariance the bispectrum will 
contain many contributions (each with a different k dependence), and it is not even well-defined what we 
mean by subleading corrections to the squeezed limit. 

Eq. (37) takes into account all the non-linearities that develop when all the modes are outside the Hubble 
radius. To this contribution we have to add the effect of the interaction among the fields which occur around 
horizon crossing, i.e. to calculate {6(p^ Scp'^ Sf^) when the modes are comfortably outside the Hubble radius. 
We want to study this correlator in the limit in which the wavevector of one of the scalars, let us call this scalar 
S(p^, is much longer than the others. Let us consider the action for all the Sip^s in spatially fiat gauge and take 
an interaction without spatial derivatives on Sip^. In the long wavelength limit we take Sip^ as homogeneous 
and the calculation reduces, similarly to the single-field case, to the calculation of the 2-point function of the 
remaining modes in the presence of this background. In this case the effect of the long mode cannot be traded 
for a coordinate rescaling, and moreover the long mode will in general be slightly time dependent because 
of its small mass; for this reason also terms with a time derivative on the long mode have to be taken into 
account^. However what is relevant is that the effect of the long mode on the short-scale 2-point function is 



''in eq. (37), we are implicitly assuming that the leading contribution to the 3-point function from eq. (37) arises 
when only one of the three legs is taken at non-linear order. It is possible, however, that this contribution vanishes and 
one has to go to higher order in Bij. For example this happens if we have two fields and (^{x) = A5(pi{x) + B5lP2{xY 
[31]; let us check that also in this case there are no linear corrections to the behaviour in the squeezed limit. The 
bispectum can be written, assuming exact scale invariance for 5(p2 and with the usual definition of kg, as 

(C(fci)C(fc2)C(fc3)> = (2^)-^<5(fci+fc2 + fc3)-6i3^A^^, / f\ \ I ^ . (38) 

J {^-^Y {q- fci/2)3 {q + ki/2Y {q + ksY 

In the limit k\ <^ ks the integral is dominated by the two divergences at if = ±A:i/2, which give a local shape with a 
logarithmic divergent coefficient [31]. To study the corrections to this, it is safe to assume g <C fcg, as the contribution 
coming from integrating outside this regime scales as ~ l/^si which is irrelevant for us. Therefore, the third power 
spectrum can be expanded in powers of q/ks as 

' '^%3i#l. (39) 



(g2-hfc|-2g-fcs)3/2 kiy 2 A;| k 

The q ■ ks term is odd if we send q — >■ —q, while the rest of the integrand is even, so we can drop it. The other terms 

just give contributions of order or higher. 

^Although the time dependence is slow-roll suppressed, we cannot neglect it as its eflects pile up out of the Hubble 

radius and can easily become of order one. 
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still proportional to the amplitude of its growing mode and therefore to the long-mode power spectrum once 
we average. This power spectrum will be ~ and, in the squeezed limit, it adds a local contribution to 

the SN calculation. 

What we now have to check is that corrections to what we said above are all quadratic in fci /ks, similarly 
to what we did in Section 2, but now in spatially flat gauge. 

1. Wavefunction corrections. As we are assuming that all the fields have a power spectrum close to scale- 
invariance, their quadratic action will be close to the one for ^, eq. (17), but with a small mass term. 
Thus the only difference will be that there is a tiny time-dependence of the mode when out of the 
Hubble radius. Corrections to this are 0{kL/ksY from the gradient terms and even more suppressed 
from the decaying mode. 

2. For terms with a spatial derivative of the long mode, the same argument as in the single- field case 
applies. 

3. Terms with a time derivative of the long mode have been considered above. Of course the time derivative 
will receive corrections from the corrections to the wavefunction, but again suppressed by at least 

0{kLlksf . 

4. Solutions for 5N and iV*. First of all, let us repeat the calculation of the single-field case in spatially 
flat gauge and then extend it to the multi-field case. The simplest way to get the constraint equations 
in the new gauge is to write the C-gauge constraints - eqs (28) and (29) - in terms of the flat gauge 
variables. Using the first order relations C = —Hit, N'^q ~ ^(flat) + diir/a^ and 5N(^q = (57V(flat) — tt", 
we get the constraints in the new gauge 

2SN{3M^iH^ + M^iH - + diN\2MlyH + ciM^) = 

(2A/|ii7 - UI^ - 3ciAI^H)Tr - ciM^V^n/a^ - 3{2M^iHH + ciM^H)tt , (40) 

{2MlyH + ciM^)di5N + (c2 + C3)A4'^d^djN^ = 

- {2Ml^H + (3c2 -f C3)HM^)d^TT + (ciA/^ - (3c2 + C3)M^H)d^7: - (cs + C3)M^d.,V^TT/a^ . (41) 

Notice that there are two small qualitative differences compared with the ^-gauge. First of all the 
solution for SN will start from a term without gradients, i.e. 0{kL/ks)'^ , plus corrections which are 
quadratic in the gradients. This does change our conclusions, but it just says that the solutions of the 
constraints now contribute to interactions terms without derivatives on the long mode. The second 
difference is that at first sight diN'^ appears of order zero in the gradients. This of course cannot occur 
because it would imply that diverges as we approach an homogeneous solution. Indeed, we know 
from the calculation in the other gauge that this does not happen and diN'' is quadratic in the gradients. 

Now let us add other fields into the game. This can be done remaining in the effective theory approach 
following [32]. In ^-gauge, additional fields S(f^ are introduced writing generic operators invariant under 
spatial but not time diffeomorphisms, for example (17'^'^ + l)g'^^df^S(f^ . Given that the 5(p^s do not enter 
in the non-linear realization of the time diffeomorphisms, it is easy to see that the new operators will 
only contribute to the right-hand side of the constraint equations in fiat gauge, with terms linear in 
the 6(p^s and their derivatives. These new terms cannot change the conclusion above that SN starts at 
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zeroth order in the gradients with quadratic corrections, while must remain linear in the gradients 
as it must go to zero as we approach an homogeneous solution. 

This concludes our proof that there are no linear corrections to the squeezed limit. We have verified this 
in the multi-DBI models of [33] and in the non-Gaussianitics generated by a cubic interaction XSkp^ of an 
isocurvature scalar [34]^. 

Notice however that we assumed from the beginning that each field has a very small mass compared to the 
Hubble scale during inflation. Indeed a clear counterexample to our conclusions is given by the compelling 
"quasi-single field" models by Chen and Wang [12]. These models feature an isocurvature direction with 
mass comparable to H and large self-coupling. This field quickly decays out of the horizon and it features a 
isocurvature power spectrum which is far from scale invariance, while the power spectrum of C, remains scale- 
invariant. However, the mixing with the adiabatic direction gives rise to a (scale-invariant) non-Gaussianity 

o . 3/2 

with a peculiar squeezed limit, ranging from kj^ , as in a local model, to fc^ , depending on the value of 
the mass. In general, once large deviations from scale-invariance are allowed, one cannot derive any general 
property for the bispectrum. 



5 The squeezing of the scale-dependent bias. 

Intuitively, since it probes the clustering of small halos on top of large-scale density perturbations, the halo 
bias is sensitive to the squeezed limit of the bispectrum [4, 5, 6], the regime studied in this paper. However, 
as we discussed in Section 3, for some models the consistency relation is valid only for a very small value of 
kh/ks- For this reason, it is important to quantify the amount of squeezing probed by realistic measurements 
of the large scale bias. In order to estimate how the halo bias depends on the bispectrum, we use the following 
expression derived in [5] for an analytical estimate of the correction to the halo bias due to primordial non- 
Gaussianity^: 



Abhik,R) 1 



bh D{z)MR{k) SttVI Jo 



dfci klMniki 



, I . F fci,^fc2 + fc? + 2fcfci/i,fc 

dA*X7v^(y'fc2 + fc2 + 2fcfci^j^ ^, (42) 

where cr^ is the variance of the dark matter density perturbations smoothed on a scale of Lagrangian radius 
R, 5c is the critical threshold for the collapse of a spherical object, which for a matter-dominated universe 
takes the value 6c = 1.686, D{z) is the linear growth factor normalized to be D{z) = (1 -I- z)~^ during matter 
domination, and F are the primordial power-spectrum and primordial bispectrum for C, and M r is the 
linear relation between the dark matter density perturbations smoothed on a scale R and the primordial 
curvature perturbation 

91.2 

-^fl(fc) ^ ^^^mH^fl(fc) , (43) 

where VLm is the present time fractional density of matter, Hq is the present Hubble rate, T is the transfer 
function normalized to one on large scales and Wr is the filter with characteristic scale R, which we will take 
here to be top-hat in real space. 

^Notice that the formula for the bispectrum given in this work contains a typo, which has been corrected in [35]. 
^Notice that in ref. [5] the factor 1/ M.R{k) is missing in their eq. (14). 
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Figure 1: Expression (42) for the modification to the halo bias in the presence of local non-Gaussianity 
multiplied by fc^ before integration over ki, for a smoothing length of i? = 2h~^ Mpc at redshift z = 1. 

Equation (42) has been shown to be in good agreement with N-body simulations [36], where several 
templates for the bispectrum were used to generate non-Gaussian initial conditions. Their results confirm 
that the scale dependent bias is sensitive to the dependence of the template on fci as it approaches the 
squeezed limit. Thus the equilateral template, which diverges like l/fci in the squeezed limit, has a negligible 
correction to the bias Abh/bh with no strong scale dependence, while the local one, which diverges like 
has a correction to the bias which grows like towards small k, and the orthogonal template, diverging 
like l/fc|_, lies in between with a correction to the bias that grows roughly like 1/k (^). This can be seen also 
from figure 1 where we plot k^Abii(k, R)/bh from equation (42) for a local bispectrum without integrating 
over ki. Notice that the plotted function goes to a constant for small k as expected, and as a function of fci it 
has a peak around fci ~ 0.2, a scale set by the transfer function for small objects {R < 10Mpc/i~^). We thus 
see that the integral in equation (42) becomes large and grows as k~^ for k much smaller than the values of 
ki contributing the most to the integral. 

Using expression (42), we can also estimate how sensitive the scale-dependent bias is to the squeezed 
limit, i.e. how small a ratio k^/ks it can probe. We do this by computing Abh{k, R) /bh for a "modified" 
template that is exactly like the local one when far from the squeezed limit, but that changes its behaviour 
as a function of k]^ from l/fc£ to I/Zcl when k]^ < akg for a given value of a. In other words, it is a local 
shape that changes to an "equilateral" shape when the ratio of the long mode to the short modes is less than 
a. At some fixed scale k and smoothing scale i?, the correction to the bias will be insensitive to this change 
in the local template for very small values of a. However, as a increases, we start deviating from the local 
case, as shown in figure 2. In order to estimate the best we can do to probe the squeezed limit, we use a scale 
of k = 0.001 /iMpc~^, which very roughly corresponds to the largest scales available in the planned EUCLID 
survey, and a smoothing length of i? = 1 Mpc /i"^, corresponding to objects with a mass of 3.1 x lO^^Mo^i"^. 
From figure 2 we see that there are 0(0. 1) changes at a « 0.001, so we expect the scale dependent bias to be 
sensitive at most to a ratio k^/ks as small as 0.001. 

*We are here referring to the standard orthogonal template defined in the main text of [37]. In the Appendix of 
the same paper another orthogonal template is defined, which goes as 1/kL in the squeezed limit. See also [20]. 
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Figure 2: Comparison of the non-Gaussian correction from the local template to the halo bias from a "modified" 
local template with the extreme squeezed limit removed (see text). The blue solid line is a "conservative" 
estimate for a smoothing scale corresponding to 2.5 x IO^^Mq/i"^, the red dashed line is an "optimistic" 
estimate for larger scales and smaller objects with a mass of 3.1 x 10^^ Mq h^^. The choice of scale for the 
"optimistic" case corresponds roughly to the largest scale accessible to the planned EUCLID survey. 

6 Conclusions 

In the last few years, a large number of models giving rise to sizeable non-Gaussianities have been proposed. 
In this context, it is useful to study general features which are present under very general assumptions. Here 
we have made more precise the single-field consistency relation for the 3-point function, showing that there 
are no linear corrections ©(fci/fcg) as we depart from the squeezed limit, once written in the form (1). The 
absence of linear corrections to the local behaviour also holds in multifield models, but only assuming that all 
the fields involved have a mass much smaller than H. These properties are a useful check for the calculation 
of the bispectrum in various models, and have obvious implications for the scale dependence of the bias. 
Single-field models do not give any scale-dependence in the bias, not surprisingly, as a long mode does not 
affect the local physics until we are sensitive to the curvature it induces, which is a 0{kj^) correction and 
corresponds to the conventional bias. Standard multifield models (containing only very light fields) can only 
give a bias going as l/kj^ (associated to the local shape), while for example it is not possible to have a bias 
going as l/k^. However, quasi-single field models can give a different scale dependence of the bias, and it 
would be interesting to study how Large Scale Structure measurements can constrain these models. 
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